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LO  lO  t> 


I. 


Introduction 


There  has  been  considerable  interest  over  the  recent  observations 
of  small-scale  irregularities  during  equatorial  Spread  F (Farley  et 
al . , 1970;  Woodman  and  LaHoz,  1976;  Costa  and  Kelley,  1978a, b;  Huba  et 
al. , 1978).  Radar  backscatter  measurements  indicate  density  fluctua- 
tions with  scale  sizes  of  3m,  lm  and  36cm.  Moreover,  since  density 
inhomogeneities  exist  with  scale  lengths  as  short  as  30m  (Costa  and 
Kelley,  1978a)  (these  scale  lengths  presumably  arise  because  of  the 
primary  growth  of  long  wavelength  fluid  modes,  e.g.,  the  collisional 
Rayleigh-Taylor  instability),  it  is  plausible  to  consider  whether  or 
not  the  excitation  of  drift  waves  could  be  responsible  for  the  small- 
scale  irregularities.  Recently,  Huba  et  al.  (1978)  have  demonstrated 
that  the  ion-drift-cyclotron  and  the  lower-hybrid-drift  instabilities 
are  viable  candidates  to  produce  the  lm  and  36cm  irregularities.  How- 
ever, the  3m  irregularities  cannot  be  explained  by  these  instabilities 
based  on  linear  theory.  Costa  and  Kelley  (1978a, b)  have  suggested  the 
universal  drift  instability  (collisionless,  low-frequency  drift  insta- 
bility) may  account  for  the  3m  irregularities.  Unfortunately,  their 
assumption  of  a collisionless  plasma  is  questionable  for 
typical  ionospheric  parameters  and  their  theory  is  limited  to  wave- 
lengths greater  than  6m. 

In  this  brief  report  we  consider  the  linear  theory  of  low 

frequency  drift  waves  for  parameters  typical  of  the  ionosphere  during 

equatorial  Spread  F.  We  use  kinetic  theory  to  describe  the  ion  and 

electron  dynamics.  (Note:  A kinetic  theory  is  required  for  the  ions 
Note:  Manuscript  submitted  February  16,  1979. 
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since  for  the  3m  irregularities  k^r^  ^ 10  >>  1,  where  r^  is  the  mean 
ion  Larmor  radius).  Electron-ion  and  electron-neutral  collisions  are 
considered  using  the  Bhatnagar-Gross-Krook  (BGK)  collision  model 
(Bhatnagar  et  al.,  1954).  We  also  include  ion- ion  collisions  via  a 
model  Fokker-Planck  equation  for  the  numerical  results  (Dougherty , 

1964;  Huba  and  Ossakow,  1979)  and  the  Landau  collision  integral  for  the 
analytic  results  (Rukhadze  and  Silin,  1969;  Mikhailovskii,  1974).  We 
find  that  a collisional  drift  wave  can  occur  in  the  3m  regime  but  that 
it  is  damped  by  ion  viscosity  unless  n < 2 x 10^  cm"^  for  density  grad- 
ient scale  lengths  ^ 75ra.  Thus,  it  is  improbable  that  low  frequency 
drift  waves  can  linearly  generate  the  3m  irregularities  observed  during 
equatorial  Spread  F. 

The  structure  of  the  paper  is  as  follows.  In  the  next  section  we 
present  the  basic  assumptions  and  the  dispersion  equation  considered 
in  the  analysis.  Section  III  contains  the  analytical  and  numerical 
results  of  our  theory.  In  the  final  section  we  discuss  the  implica- 
tions of  our  results  concerning  the  3m  irregularities. 

II . Basic  Assumptions  and  the  Dispersion  Equation 

The  geometry  we  consider  is  described  as  follows.  The  ambient 
magnetic  field  is  JS  = BQez,  the  density  depends  only  on  the  x coordi- 
nate (n  = n^(x))  and  the  temperature  is  assumed  to  be  constant.  Each 

+ 

species  a(electrons  and  0 ions)  has  a drift  velocity  V = e^ 

where  = (v  2/2ft^)  d Jln/dx  is  the  diamagnetic  drift  velocity, 

v = (2T  /m  )**  is  the  thermal  velocity,  Q = e B /m  c is  the  cyclotron 
a a or  J a a o a 

frequency  and  n = n^  ~ n_^.  A net  current  exists  in  the  plasma 
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J = en(V,  -V , )e  = enV,,(l  4-  T /T,)e  which  provides  the  free 
nuO  v di  de  y di  e i y r 

energy  to  drive  an  instability.  Perturbed  quantities  vary  as 
exp  i(k^y  4-  k z - cot)  and  we  consider  only  electrostatic  oscilla- 
tions. We  make  use  of  the  local  approximation  which  requires  kL^  >>  1 
where  = (d  In  n/dx)-1  is  the  density  gradient  scale  length.  We 

include  electron-ion  and  electron-neutral  collisions  (v  = v , 4-  v ) 

e ei  en7 

and  ion-ion  collisions  - v ).  We  ne8lect  ion-neutral  collisions 

by  requiring  « kj v^  which  restricts  the  results  of  the  analysis  to 

2 2 

altitudes  above  300  km.  Finally,  we  consider  u>  <<  ft  , Q , k r_  , >>  1 

i e Li 

2 2 

k rT  <<  1 and  v , , <<  . 

Le  ii  i 

The  dispersion  equation  is  found  to  take  the  following  form  within 
the  context  of  our  model 

0 (1) 


D(w,k)  = 1 + xe+  X± 


where  (Kadomtsev,  1965) 


2(J  2 

X = 

e k2v  2 


1 -k  V , . 

1 + v x de  Z(C  ) { 

k||Ve  e * 


1 + 


iv 

k^-  Z^e> 


— 1 

} ' 


(2) 


and  (Huba  and  Ossakow,  1979) 


2i'  2 f 

X±  = —£L~ 

kV  L 


io -k  V..  ® 

1 + i / dr  exp  (i  t-^Ct) 

i o i 


(3) 


Here , 


[ i ~V,T  "I  kip  b r -v,t-i 

cos  0 4-  v,x-e  cos(a-0)  4-  v,  t-1  4-  e 

1 J v L 1 J 


k V,,  (■  -V . T -V . T -j 

+ i 1 le  sinr  + v^l-e  cost)I  , 
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£e  = (<jo+ive)/k|jve,  a)pa2  = 4ime2/ma,  v±  = 6 = 2tan  1 v±  and 

t>i  = ^k±2r^2.  Several  comments  are  in  order  concerning  the  validity 
of  this  dispersion  equation.  Equation  (2)  is  based  on  the  BGK 
collision  model  and,  strictly  speaking,  only  describes  electron-neutral 
collisions.  However,  this  model  reproduces  the  results  based  on  fluid 
theory  in  the  collisional  regime  and  the  kinetic  results  in  the  colli- 
sionless regime.  Here  we  define  collisional  (collisionless)  as 

>>  k*  ve(ve  <<:  k ve)  • Thus,  only  in  the  semi-collisional  regime 
(v^  ~ kjjv^)  is  Eq.  (2)  suspect  and  this  regime  is  not  relevant  for  the 
ionospheric  parameters  of  interest.  Equation  (3)  is  based  upon  a 
model  Fokker-Planck  equation  (Dougherty,  1964)  and  is  valid  only  for 
Vii  <<:  ^i  anc^  kj_2rLi2  >:>  1 which  is  precisely  the  regime  of  interest. 
The  numerical  results  presented  in  the  next  section  will  be  based  on 
Eqs.  (1)- (3) . 

Although  Eq.  (3)  does  not  lend  itself  to  a simpler  analytical 

expression  for  the  parameters  we  are  concerned  with,  an  alternative 

formulation  of  X—  based  upon  the  Landau  collisional  integral  (Rudkadze 

and  Silin,  1969)  permits  analytic  solutions  to  Eq.  (1)  in  several 

limits.  In  particular,  for  k 2r  , 2 >>  1,  u>  » ki  v.  and 

± Ll  ||  'i 

v^k±2r^2  <<  w one  can  show  that  (Mikhailovskii , 1974) 


X 

i 


2w2 


k2v^2 


1 - 


k V, . 


/IT 


kxrLi 


(1-i 


3(ir+l) 

4/2 


-^k2r  2) 
ui  x Li' 


(4) 


We  will  use  this  expression  in  the  following  section  to  derive  simple 
expressions  for  the  eigenf requency , w = or  + iy. 
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III.  Theoretical  Analysis 


A.  Analytical  results. 

We  now  present  several  analytical  results  for  a variety  of 
conditions  to  elucidate  the  physical  nature  of  the  instability  and  the 
effects  of  collisions.  For  all  the  cases  considered  we  will  assume 


that  k±2rLi2  » 1,  k||Ve  » o>  » k|jvi,  to  « k^,  T0  = T±, 

2v  2 


k ^rT  * <<  1,  <<  Q , , Q and  y <<  m . 

± Le  i’  e 1 r 


1.  v..  = 0 and  v =0 

ii  e 

In  the  collisionless  limit  the  dispersion  equation  becomes 

k V. 


D(aj,k)  « 2- 


k V, 

-L  de  1 


r~  k r 
/tt  ± Li 


- i /t\ 


± de 
k,i  v 


= 0 


(5) 


which  yields  the  solution 


k V . 

± de 


0)  = 

r 


2/tt  k rT  , 
± Li 


Y ■ 


k 2V,  2 
± de 

2kiiv 


k rT  , 
± Li 


6) 


We  comment  that  the  instability  is  kinetic  and  is  excited  through  an 
electron-wave  resonance.  Relaxing  the  restriction  y <<  it  is  found 


that  (Mikhailovskii,  1974) 

..  i v- 

Y\ 


'M 


to  « 

r 8/n 


(7) 


for  kn  « tt- 


/tt 

2 


k V, 

J.  de 


where  y^  denotes  the  maximum  growth  rate. 


It  is 


. e 

interesting  to  note  that  in  the  limit  k^2r^2  >>  1,  y^  is  independent 
of  wavenumber. 
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2.  v..  ^ 0 and  v = 0 

ii  e 


In  this  limit,  in  which  we  neglect  electron  collisions  but 
include  ion- ion  collisions,  the  dispersion  equation  becomes 

k'v‘-  ’ ■ ' 1 1WI  vli  t 


D(co,k)  = 2 


k V , , , / 

-i  de  1 1 ( i _ 

u)  r~  k rT  . \ 

/it  i Li  ' 


4/2 


c 2r  2) 

± Li  ) 


k va 

- 1 0 

kl|V 


(8) 


where  we  have  assumed  v^kA2r^2  <<:  Solving  Eq.  (8)  yields 


03 

r 


k.V, 

-l  de 

2/iTk  rT  . 
-l  Li 


k 2V,  2 
± de 

2kii  v 


kxrLi 


3 ( TT+1 ) 

4 Jl 


v k 2r  “ 
Vlik±  rLi 


(9) 


Again,  the  instability  is  kinetic  but  ion  viscosity  provides  a damping 

mechanism  which  is  particularly  effective  at  the  shorter  wavelengths. 

3.  v,.  i~  0 and  v f 0 

ii  e 

We  finally  consider  the  fully  collisional  limit  and  assume 

v..k  2rT  2 <<  o)  03  <<  v , v >>  kiiy  and  v o)  < kn2v  2.  The  disper- 
li  ± Li  e e ||  e e ||  e r 

sion  equation  is 


D(03,k)  = 2-^^  — ~ — (l  “ 

Cl)  r-  k rT  , \ 

/tt  i Li  \ 


i 3 (tt+1)  Vii  k 
4/2  W 


2r  2\ 
' ± Li  ) 


- 2i 


v k V. 
e i de 

k||2ve2 


(10) 


for  which 


0) 

r 


ki.Vde 

2^kxrLi 


; y = 


VekxVde  _ 3 (tt+1)  2 2 

4k  2ve2  4/2  11  A L1  ' 


(11) 


6 


In  this  limit  the  nature  of  the  instability  has  changed  from  being  a 
kinetic  (collisionless)  instability  to  a collisional  instability. 
However,  ion  viscosity  is  still  an  effective  damping  mechanism.  For 
parameters  typical  of  the  ionosphere  this  final  limit  is  the  most 
relevant.  It  should  be  noted  that  many  of  the  approximations  made 
which  led  to  the  simplified  dispersion  equations  breakdown  for 
realistic  values  of  v^,  and  and  a numerical  analysis  is  required 
to  determine  w. 

B.  Numerical  results. 

We  now  present  the  results  of  a numerical  analysis  based  upon 

Eqs.  (l)-(3).  The  parameters  chosen  which  are  common  to  both  Figs.  (1) 

and  (2)  are:  T /T,  = 1.0,  U)  /P  = 10.0  and  V,,/v,  = 0 . 0A  which 
e i pe  e di  i 

corresponds  to  a density  gradient  scale  length  L ~ 75ra  for  a 1000°K 

n 

0+  plasma  in  a 0.3  Gauss  field.  We  mention  the  results  are  relatively 
insensitive  to  the  parameter 

In  Fig.  1 we  present  the  growth  rate  spectrum  y/P^  vs.  k^r^ ^ 

for  different  collisional  conditions.  The  solid  curve  ( — ) represents 

the  collisionless  case  (v^  = 0 and  = 0) , the  dashed  curve  ( ) 

represents  only  electron  collisions  (v^  = 0 and  = 2.0),  the 

dotted  curve  ( ) represents  only  ion-ion  collisions  - 1-0 

x 10  ^ and  = 0)  and  the  dash-dot  curve  (-•-*-•)  is  fully  collisional 

(Vn/Pi  = 1.0  x 10  5 and  v^/P^  * 2.0).  The  value  of  kj/k^  is  chosen  to 

give  maximum  growth  for  each  value  of  k r . Also,  as  mentioned 

i-  Li 

earlier,  Eq.  (3)  is  only  valid  for  k ?r  . 2 >>  1 when  v ^ 0 and  has 

* ^ Ll  11 

been  used  for  k^r^  > 4.0.  For  k^r^  < 4.0  the  role  of  ion  viscosity 
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becomes  negligible  and  the  spectra  resemble  the  v^“Q  spectra  in  this 
region.  Several  points  of  interest  concerning  these  curves  are  the 
following.  First,  maximum  growth  for  the  collisionless  case  occurs,  for 
k,rTJ  w 0.8  which  is  in  agreement  with  the  results  of  Gary  and 

L»1 

Sanderson  (1978).  We  mention  that  Costa  and  Kelley  (1978b)  incorrectly 
found  maximum  growth  for  k^r^  w 1.5.  Electron  collisions  cause  the 
peak  growth  to  shift  to  higher  k (kxrn  w 1.0)  while  ion-ion 
collisions  do  not  affect  the  wavenumber  of  maximum  growth.  Second, 
for  k^r^ £ > 4 the  growth  rate  is  relatively  constant  for  a large  range 
of  as  anticipated  from  Eq.  (7)  for  the  collisionless  case. 

Third,  the  electron  collisional  drift  instability  has  lower  growth 
rates  over  most  of  the  spectrum  (except  for  k 2rT  ,2<<  1)  than  the 
collisionless  instability.  Finally,  ion-ion  collisions  heavily  damp 
the  short  wavelength  modes  and  can  stabilize  the  instability. 

In  Fig.  2 we  plot  y /ft,  vs.  v /ft.  for  v,,/ft,  * 0.,  1.0  x 10  ^ , 

ml  el  ll  l 

2.0  x 10~5,  3.0  x 10-5,  4.0  x 10"5,  5.0  x 10-5,  and  k rL±  = 10.0 
(which  roughly  corresponds  to  3m  wavelengths).  Here  y^  denotes  the 
maximum  value  of  the  growth  rate  with  respect  to  kji/k^.  The  important 
features  of  this  curve  are:  (1)  electron  collisions,  although  res- 
ponsible for  the  instability  when  vg/ft^  > 0.6,  have  a stabilizing 
influence  as  is  increased  and  (2)  relatively  infrequent  ion-ion 
collisions  can  stabilize  the  instability  due  to  ion  viscuous  damping. 

IV.  Discussion 

We  have  presented  an  analysis  of  low  frequency  drift  waves  based 
on  kinetic  theory  including  electron  (i.e.,  electron-ion  and  electron- 
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neutral)  and  ion  (i.e.,  ion-ion)  collisions.  The  primary  purpose  is 
to  determine  whether  or  not  drift  waves  can  produce  the  3m  irregulari- 


ties observed  during  equatorial  Spread  F based  on  linear  theory.  We 
now  address  this  question. 

For  typical  ionospheric  parameters  it  is  expected  that 

v ~ 200-1000  sec  ^ and  ft  200  sec  ^ so  that  v /ft,  ~ 1-5.  Based  on 
e i e i 

our  theory  we  find  that  a collisional  drift  wave  can  be  excited  in 

this  parameter  range.  However,  ion  viscosity  can  heavily  damp  the 

short  wavelength  modes  (k  ?rT  2 >>  1)  and  even  infrequent  ion-ion 

collisions  can  stabilize  the  instability.  From  Fig.  2 we  see  that  the 

3m  mode  is  stabilized  for  v,,/ft,  > 5 x 10  ^ when  L ~ 75m  and  v > ftJ . 

ii  i n e i 

For  = 1000°K  and  B = 0.3  Gauss  this  corresponds  to  n > 2 x 10^  cm  ^ 
which  is  generally  the  situation  in  the  F region.  The  exception,  of 
course,  would  be  deep  inside  a plasma  depletion  (bubble)  where  plasma 
densities  can  fall  below  this  value.  Still,  this  would  not  explain  the 
3m  backscatter  return  observed  on  the  bottomside  of  the  F region  by 
Woodman  and  LaHoz  (1976)  just  after  sunset  when  Spread  F begins. 

An  approximate  analytic  expression  can  be  derived  relating  density 
and  density  gradient  scale  length  to  determine  the  conditions  under 
which  instability  can  occur.  We  demand  that  > v^k^2rj  ^2/2  for 
instability  where 


Yo  = L 


LLi 


8/tt 


(12) 


is  the  maximum  growth  of  the  collisionless  instability  in  the  regime 

k 2r  2 >>  i (Eq.(7)).  Although  the  collisionless  mode  is  usually  not 
X Ll 

excited,  its  growth  rate  is  not  significantly  different  from  the 
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collisional  mode.  For  = 1000°K  and  = 200  sec-^  we  obtain  the 
following  condition  for  instability 


„ Li  1 „ ,„6  -3 

t~  — 7 2 x 10  cm  • 


(13) 


3 -3 

Thus,  for  k r_  . = 10  and  L /rT 4 ^ 10  we  require  n < 2 x 10  cm  for 
i Li  n Li  M 


instability  which  is  in  good  agreement  with  the  numerical  estimates. 

In  conclusion,  it  seems  unlikely  that  a linear  theory  of  drift 

waves  can  account  for  the  observation  of  3m  irregularities  during 

equatorial  Spread  F.  Although  it  is  possible  to  generate  sharp  density 

gradients  (L  ~ 10  rT  . ) on  the  edges  of  density  depletions  (bubbles) 
n li 

which  are  sufficient  to  excite  a collisional  drift  instability,  the 
condition  (Eq.  (13))  to  overcome  ion  viscuous  damping  is  difficult.  On 
the  other  hand,  longer  wavelength  modes  (k  rT  . . ^ 1)  are  more 
readily  excited  since  ion  viscuous  damping  is  less  effective  for  these 
modes.  We  suggest  it  may  be  possible  for  a large  amplitude,  long 
wavelength  mode  to  nonlinearly  generate  the  short  wavelength  turbulence 
via  a parametric  process  and  we  are  presently  investigating  this 
possibility. 


Ac  knowl ed  g emen  t s 


We  would  like  to  thank  Dr.  P.  K.  Chaturvedi  for  several  helpful 
discussions.  This  research  has  been  supported  by  the  Defense  Nuclear 
Agency  and  the  Office  of  Naval  Research. 


10 


References 


Bhatnagar,  P.  L. , E.  P.  Gross  and  M.  Krook,  A Model  for  Collision 
Processes  in  Gases.  I.  Small  Amplitude  Processes  in 
Charged  and  Neutral  One-Component  Systems,  Phys.  Rev.,  94 , 
511,  1954. 

Costa,  E.  and  M.  C.  Kelley,  On  the  Role  of  Steepened  Structures 

and  Drift  Waves  in  Equatorial  Spread  F,  J,  Geophys.  Res.,  83 , 
4359,  1978a. 

Costa,  E.  and  M.  C.  Kelley,  Linear  Theory  for  the  Collisionless 
Drift  Wave  Instability  with  Wavelengths  Near  the  Ion 
Gyroradius,  J.  Geophys.  Res.,  83,  4365,  1978b. 

Dougherty,  J.  P.,  Model  Fokker-Planck  Equation  for  a Plasma  and 
Its  Solution,  Phys.  Fluids,  1,  1788,  1964. 

Farley,  D.  T.,  B.  B.  Balsley,  R.  F.  Woodman  and  J.  P.  McClure, 

Equatorial  Spread  F:  Implications  of  VHF  Radar  Observations, 
J.  Geophys.  Res.,  75 , 7199,  1970. 

Gary,  S.  P.  and  J.  J.  Sanderson,  Density  Gradient  Drift 

Instabilities:  Oblique  Propagation  at  Zero  Beta,  Phys . 

Fluids,  21,  1181,  1978. 

Huba,  J.  D.,  P.  K.  Chaturvedi,  S.  L.  Ossakow  and  D.  M.  Towle,  High 
Frequency  Drift  Waves  with  Wavelengths  Below  the  Ion  Gyro- 
radius  in  Equatorial  Spread  F,  Geophys.  Res.  Lett.,  5^  695, 
1978. 

Huba,  J.  D.  and  S.  L.  Ossakow,  Destruction  of  Cyclotron 

Resonances  in  Weakly  Collisional,  Inhomogeneous  Plasmas,  (to 
be  published  in  Phys.  Fluids,  1979) , 


11 


Kadomtsev,  B.  B. , Plasma  Turbulence,  Academic  Press,  N.  Y.,  1965 


Mikhailovskii,  A.  B.,  Theory  of  Plasma  Instabilities  Vol.  II, 
Consultants  Bureau,  N.  Y.,  1974. 

Rukhadze,  A.  A.  and  V.  P.  Silin,  Kinetic  Theory  of  Drift- 

Dissipative  Instabilities  of  a Plasma,  Sov.  Phys . Uspekhi, 
11,  659,  1969. 

Woodman,  R.  F.  and  C.  LaHoz,  Radar  Observations  of  F-Region 

Equatorial  Irregularities,  J.  Geophys.  Res.,  81,  5447,  1976 


12 


y/ft,  ( X 10"’) 


Figure  1.  Growth  rate  spectrum  v/ft.  vs.  k rT  . for  T / T.  = 1.0, 

i x Li  e i 

I 

ajpe/^e  = 10*0  and  = 0,04  (L^  « 75m).  The  growth 

rate  has  been  maximized  with  respect  to  k for  each  value 

-4  -3 

of  and  typically  /k^  ^10  - 10  Here,  the 

collisional  parameters  are: 


( — 

II 

1 

0;  v 

e 

= 0 

( — 

— ) V = 

0;  v 

= 2.0  0, 

ii 

e 

i 

<•••■ 

vi±  = 

1 . 0 x 

10-5  n.;  v 
i e 

II 

•H 

1 

1.0  x 

io-5  n.,  v 

i e 

13 


01  X)  !u/ nx 


6.0 


Figure  2.  A plot  of  y /ft.  vs.  v /ft.  for  the  same  parameters  as 

m 1 e i 

Fig.  1 with  = 10.0  and  = 1*0  x 10  \ 
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NATIONAL  TELECOMMUNICATIONS  & INFO  ADMIN 
BOULDER,  CO  80303 

01CY  ATTN  A.  JEAN  (UNCLASS  ONLY) 
01CY  ATTN  W.  UTLAUT 
0 ICY  ATTN  O.  CROMBIE 
0 1 CY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  & ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTS  NT  OF  COMMERCE 
BOULDER,  CO  80302 

01CY  ATTN  R.  GRUBB 

01CY  ATTN  AERONOMY  LAB  G.  REIO 


OEPARTMENT  OF  DEFENSE  CONTRACTORS 


AEROSPACE  CORPORATION 

P.  O.  BOX  92957 

LOS  ANGELES,  CA  90009 

01 CY  ATTN  1.  GARFUNKEL 
0 ICY  ATTN  T.  SALMI 
01CY  ATTN  V.  JOSEPHSON 
01CY  ATTN  S.  BOWER 
01CY  ATTN  N.  STOCKWELL 
0 ICY  ATTN  0.  OLSEN 
0 ICY  ATTN  U.  CARTER 
0 ICY  ATTN  F.  MORSE 
01CY  ATTN  SMFA  FOR  PWW 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5 OLO  CONCORD  ROAD 
BURLINGTON,  MA  01803 

0 ICY  ATTN  RAOIO  SCIENCES 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.  0.  BOX  983 
BERKELEY,  CA  94701 

0 ICY  ATTN  J.  WORKMAN 


BOEING  COM3 ANY,  THE 
P.  0.  BOX  3707 
SEATTLE,  WA  98124 

01CY  ATTN  G.  KEISTER 
01CY  ATTN  D.  MURRAY 
01CY  ATTN  G.  HALL 
01CY  ATTN  J.  KEM^EY 

CALIFORNIA  AT  SAN  DIEGO,  UNIV  OF 
IPAPS,  B-019 
LA  JOLLA,  CA  92093 

01CY  ATTN  HENRY  G.  BOCXER 

3RCWN  ENGINEERING  COMPANY,  INC. 

CUMMI NGS  RESEARCH  PARK 
HUNTSVILLE,  Al  35807 

01CY  ATTN  ROMEO  A.  DELIBERIS 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 

555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 

01CY  ATTN  0.  B.  COX 
01CY  ATTN  J.  P.  GILMORE 

COMPUTER  SCIENCES  CORPORATION 
6565  ARLINGTON  BLVD 
FALLS  CHURCH,  VA  22046 
01CY  ATTN  H . BLANK 
01CY  ATTN  JOHN  SPOOR 
01CY  ATTN  C.  NAIL 

COMSAT  LABORATORIES 
Ll NTH  I CUM  ROAD 
CLARKSBURG,  MD  20734 

01CY  ATTN  G.  HYOE 

CORNELL  UNIVERSITY 

OEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

01CY  ATTN  D.  T.  FARLEY  JR 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

R I CHAR  OS  ON,  TX  75080 

OICY  ATTN  H.  LOGSTON 

01CY  ATTN  SECURITY  (PAUL  PHILLIPS) 

ESL  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94066 

OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 
OICY  ATTN  C.  W.  PRETT1E 

FORO  AEROSPACE  S COMMUNICATIONS  CORP 
3939  FABIAN  WAY 
PALO  ALTO,  CA  94303 

OICY  ATTN  J.  T.  MATTINGLEY 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GOOOARD  BLVD  KING  OF  PRUSSIA 
P.  0.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.  H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.  0.  BOX  1122 
SYRACUSE,  NY  13201 

OICY  ATTN  F.  REIBERT 

GENERAL  ELECTRIC  COMPANY 
TEYPO-CENTER  FOR  AOVANCEO  STUDIES 
816  STATE  STREET  (P.O,  ORAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OICY  ATTN  DASIAC 
OICY  ATTN  DON  CHANDLER 
OICY  ATTN  TOM  BARRETT 
OICY  ATTN  TIM  STEPHANS 
OICY  ATTN  wARREN  S.  KNAPP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GAMBILL 
OICY  ATTN  MACK  STANTON 
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GENERAL  electric  tECH  SERVICES  CO.,  InC. 
nMES 

COURT  STREET 
SYRACUSE,  NY  13201 

01CY  ATTN  G.  MIlLMAN 

GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  OIVISION 
P.  0.  BOX  6770 
SANTA  BARBARA,  CA  93111 

01CY  ATTN  JOHN  ISE  JR 
0 ICY  ATTN  JOEL  GARBARINO 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

Call  class  attn:  security  officer) 

01CY  ATTN  T.  N.  OAVIS  CUNCL  ONLY) 

0 ICY  ATTN  NEAL  BROWN  (UNCL  ONLY) 

0 ICY  ATTN  TECHNICAL  LIBRARY 

GTE  SYL VANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  OIV 
77  A STREET 
NEEDHAM,  MA  02194 

0 ICY  ATTN  MARSHAL  CROSS 

ILLINOIS,  UNIVERSITY  OF 

OE  PART  ME  NT  OF  ELECTRICAL  ENG  I BERING 

URBANA,  IL  61803 

0 ICY  ATTN  K.  YEH 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
801  S.  WRIGHT  STREET 
URBANA,  IL  60680 

(ALL  CORRES  ATTN  SECURITY  SUPERVISOR  FOR) 
01CY  ATTN  K.  YEH 

INSTITUTE  FOR  OEFENSE  ANALYSES 
400  ARMY  -NAVY  ORIVE 
ARLINGTON,  VA  22202 

01CY  ATTN  J.  M.  AE IN 
0 ICY  ATTN  ERNEST  BAUER 
0 ICY  ATTN  HANS  WOLFHARO 
0 ICY  ATTN  JOEL  BENGSTON 

HSS,  INC. 

2 ALFREO  CIRCLE 
BEDFORD,  MA  01730 

01CY  ATTN  DONALO  HANSEN 

I NTL  TEL  C TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

0 ICY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

1401  C AMI  NO  DEL  MAR 
OEL  MAR,  CA  92014 

01CY  ATTN  S.  R.  GOLDMAN 


JOHNS  HOPKINS  UNIVERSITY 
APPLIEO  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

01CY  ATTN  DOCUMENT  LIBRARIAN 
01CY  ATTN  THOMAS  POTEM*A 
01CY  ATTN  JOHN  OASSOULAS 

LOCKHEEO  MISSILES  l SPACE  CO  INC 
P.  O.  BOX  504 
SUNNYVALE,  CA  94088 

01CY  ATTN  OEPT  60-12 
01CY  ATTN  D.  R.  CHURCHILL 

LOCKHEED  MISSILES  AND  SPACE  CO  INC 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

01CY  ATTN  MARTIN  WALT  DEPT  52-10 
0 ICY  ATTN  RICHARO  G.  JOHNSON  OEPT  52-12 
01CY  ATTN  W.  L.  I M"tOF  OEPT  52-12 


KAMAN  SCIENCES  CORP 
P.  O.  BOX  7463 
COLORADO  SPRINGS,  CO  80933 
01CY  ATTN  T.  MEAGHER 

L1MCAB1T  CORP 
10453  ROSELLE 
SAN  OIEGO,  CA  92121 

01CY  ATTN  IRWIN  JACOBS 

LOVELL  RSCH  FOUNDATION,  UNIVERSITY  OF 
450  AIKEN  STREET 
LOWELL,  Mft  01854 

01CY  ATTN  K.  BI8L 

M.I.T.  LINCOLN  LABORATORY 
P.  O.  BOX  73 
LEXINGTON,  MA  02173 

OlCY  ATTN  DAVIO  M.  TOWLE 
01CY  ATTN  P.  WALDRON 
OlCY  ATTN  L.  LOUGHLIN 
OlCY  ATTN  0.  CLARK 

MARTIN  MARIETTA  CORP 
ORLANDO  OIVISION 
P.  0.  BOX  5837 
ORLANDO,  FL  32805 

OlCY  ATTN  R.  HEFFNER 

MCDOWELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENLE 
HUNTINGTON  BEACH,  CA  92647 
OlCY  ATTN  N.  HARRIS 
OlCY  ATTN  J.  MOULE 
OlCY  ATTN  GEORGE  MR 02 
OlCY  ATTN  W.  OLSON 
OlCY  ATTN  R.  W.  HALPRIN 
OlCY  ATTN  TECHNICAL  LIBRARY  SERVICES 


MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  93101 


OlCY 

ATTN 

P.  FISCHER 

OlCY 

ATTN 

w.  F.  CRtVIER 

OlCY 

ATTN 

STEVEN  L.  GUTSCHE 

OlCY 

ATTN 

D.  SAPPENFIELO 

OlCY 

ATTN 

R.  BOGUSCH 

OlCY 

ATTN 

R.  HENDRICK 

OlCY 

ATTN 

RALPH  KILB 

OlCY 

ATTN 

OAVE  SCVLE 

OlCY 

ATTN 

F.  FAJEN 

OlCY 

ATTN 

M.  SCHE1BE 

OlCY 

ATTN 

CONRAD  L.  LONGMIRE 

OlCY 

ATTN 

WARREN  A.  SCHLUETER 

MITRE  CORPORATION,  THE 

P.  0.  BOX  208 

BEDFORO,  MA  01730 

OlCY 

ATTN 

JOHN  MORGANS  TERN 

OlCY 

ATTN 

G.  HAROING 

OlCY 

ATTN 

C.  E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
18  20  DOLLY  MAO I SON  BLVO 
MCLEAN,  VA  22101 

OlCY  ATTN  W.  HALL 
OlCY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVE RF I ELO  BLVO. 

SANTA  MONICA,  CA  90404 

OlCY  ATTN  E.  C.  F I ELO  JR 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 

(NO  CLASS  IF IEO  TO  THIS  ADORESS) 

OlCY  ATTN  IONOSPHERIC  RESEARCH  LAB 
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PHOTOMETR I CS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

01CY  ATTN  IRVING  L.  K OF SKY 

PHYSICAL  DYNAMICS  INC. 

P.  D.  8DX  3027 
BELLEVUE,  WA  98009 

01CY  ATTN  E.  J.  FREMOUW 


SCIENCE  APPLICATIONS,  INCORPORATED 
8400  WE ST PARK  DRIVE 
MCLEAN,  VA  22101 

01 CY  ATTN  J.  COCKAYNE 

SCIENCE  APPLICATIONS,  INC. 

80  MISSION  DRIVE 
PLEASANTON,  CA  94566 
01 CY  ATTN  S 2 


PHYSICAL  DYNAMICS  INC. 

P.  D.  BOX  1069 
BERKELEY,  CA  94701 

01 CY  ATTN  A.  THOMPSON 

R ( D ASSOCIATES 

P.  D.  BOX  9695 

MARINA  DEL  REY,  CA  90291 

01CY  ATTN  FORREST  GILMORE 
0 1 CY  ATTN  BRYAN  GABBARD 
01CY  ATTN  WILLIAM  8.  WRIGHT  UR 
01 CY  ATTN  ROBERT  F.  LELEV1ER 
0ICY  ATTN  WILLIAM  U.  KARZAS 
01CY  ATTN  H.  DRY 
01CY  ATTN  C.  MACDONALD 
01CY  ATTN  R.  TURCO 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

0 1 CY  ATTN  CULLEN  CRAIN 
01CY  ATTN  ED  8EDROZ1AN 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

01CY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 

P.  D.  BOX  2351 
LA  JOLLA,  CA  92038 

01CY  ATTN  LEWIS  M.  LlNSON 
01CY  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  D.  SACHS 
01CY  ATTN  E.  A.  STRAKER 
0 ICY  ATTN  CURTIS  A.  SMITH 
01CY  ATTN  JACK  MCDOUGALL 

RAYTHEON  CD. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01CY  ATTN  BARBARA  ADAMS 


SRI  INTERNATIONAL 
333  RAVENS WOOO  AVENUE 
MEN-D  PARK,  CA  94025 

01 CY  ATTN  DONALD  NEILSON 
01CY  ATTN  ALAN  BURNS 
0 ICY  ATTN  G.  SMITH 
01CY  ATTN  L.  L.  CDBB 
01CY  ATTN  DAVID  A.  JOHNSON 
0ICY  ATTN  WALTER  G.  CHESNUT 
01CY  ATTN  CHARLES  L.  RINO 
01 CY  ATTN  WALTER  JAYE 
01 CY  ATTN  M.  BARON 
01 CY  ATTN  RAY  L.  LEADABRAND 
01CY  ATTN  G.  CARPENTER 
01CY  ATTN  G.  PRICE 
01CY  ATTN  J.  PETERSON 
01 CY  ATTN  R.  HAKE,  JR. 

01 CY  ATTN  V.  GONZALES 
01CY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CDRP 
75  WIGGINS  AVENUE 
BEOFORD,  MA  01730 

01CY  ATTN  W.  P.  BOQUIST 

TRW  DEFENSE  £ SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONDO  BEACH,  CA  90278 

0ICY  ATTN  R.  K.  PLEBUCH 
01CY  ATTN  S.  ALTSCHULER 
01CY  ATTN  D.  DEE 

VISIDYNE,  INC. 

19  THIRD  AVENUE 

NORTH  WEST  INOUSTRIAi.  PARK 

BURLINGTON,  MA  01803 

01CY  ATTN  CHARLES  HU^HREY 
01CY  ATTN  J.  W.  CARPENTER 


SCIENCE  APPLICATIONS,  INC. 
HUNTSVILLE  DIVISION 
2109  W.  CLINTON  AVENUE 
SUITE  700 

HUNTSVILLE,  AL  35805 

01CY  ATTN  DALE  H.  DIVIS 
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IONOSPHERIC  MOOELING  DISTRIBUTION  UST 
UNCLASS  IF IEO  ONLY 


PLEASE  01  STRI  suite  ONE  COPY 


TO  EACH  OF  THE  FOLLOWING  PEOPLE: 


ADVANCEO  RESEARCH  PROJECTS  AGENCY  (ARPA) 
STRATEGIC  TECHNOLOGY  OFFICE 
ARLINGTON,  VIRGINIA 


HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MASS.  02138 


CAPT.  DONALO  M.  LEVINE 


DR.  M.  B.  MCELROY 
DR.  R.  L IND2EN 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 

OR.  P.  MANGE 
DR.  R.  MEIER 

DR.  E.  S2US2C2EWIC2  - COOE  7127 


OR.  J.  GOODMAN  - CODE  7560 


SCIENCE  APPLICATIONS,  INC. 
1250  PROSPECT  PLA2A 
LA  JOLLA,  CALIFORNIA  92037 

DR.  O.  A.  HAMLIN 
OR.  L.  LINSON 
DR.  O.  SACHS 


DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL  LABORATORY 
NOAA 

BOULDER,  COLORADO  80302 

OR.  A.  GLEW  JEAN 
DR.  G.  w.  ADAMS 
OR.  D.  N.  ANDERSON 
OR.  K.  OAVIES 
DR.  R.  F.  DOWELL Y 


A.  F.  GEOPHYSICS  LABORATORY 
L.  G.  HANSOM  FIELD 
BEDFORD,  MASS.  01730 

OR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R . SAGALYN 
DR.  J.  M.  FORBES 
DR.  T.  J.  KENESHEA 
DR.  J.  AARONS 


OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VIRGINIA  22217 

DR.  H.  MULLANEY 


PEWSTL VANIA  STATE  UNIVERSITY 
JNIVERSITY  PARK,  PEWSYLVANIA  16802 


DR. 

J. 

S.  NISBET 

OR. 

P. 

R.  RCHRBAUGH 

OR. 

D. 

E.  BAPAN 

DR. 

L. 

A.  CARPENTER 

DR.  M. 

LEE 

OR. 

R. 

OlVANY 

DR.  P. 

BENNETT 

DR. 

E. 

klevans 

UNIVERSITY  OF  CALIFORNIA,  LOS  ANGELES 

405  HILLGARO  AVENUE 

LOS  ANGELES,  CALIFORNIA  90024 

DR.  F.  V.  COR ON  I T I 
OR.  C.  KENNEL 


UNIVERSITY  OF  CALIFORNIA,  BERKELEY 
BERKELEY,  CALIFORNIA  94720 

DR.  M.  HUDSON 


UTAH  STATE  UNIVERSITY 
4TH  N.  AND  8TH  STREETS 
LOGAN,  UT*i  84322 

OR.  P.  M.  0ANCS 
DR.  R.  HARRIS 
OR.  V.  PETERSON 
OR.  R.  MEGILL 
DR.  K.  BAKER 

CORNELL  UNIVERSITY 
ITHACA,  NEW  YORK  14850 

DR.  W.  E.  SwART 2 
OR.  R.  SUOAN 
DR.  O.  FARLEY 
DR.  M.  KELLEY 
DR.  E.  OTT 


NASA 

GOOCARO  SPACE  FLIGHT  CENTER 
GREEN8ELT,  MARYLAND  20771 


COMMANDER  DR.  S.  CHANORA 

NAVAL  ELECTRONICS  LABORATORY  CENTER  OR.  K.  MAEDO 

SAN  DIEGO,  CALIFORNIA  92152 


DR.  M.  BLEIWEISS 
DR.  I.  ROTI^ULLER 
DR.  V.  HILDEBRAND 
MR.  R.  ROSE 


PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  NEW  JERSEY  08540 

OR.  F.  PERKINS 
DR.  E.  FRIEMAN 


U.  S.  ARMY  ABERDEEN  RESEARCH  AND  DEVELOPMENT  CENTER 


BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MARYLANO 

DR.  J.  HE  I ME  RL 


COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  O.C.  20360 

DR.  T.  C2U8A 


INSTITUTE  FOR  OEFcNSE  ANALYSIS 
400  ARMY /NAVY  ORIVE 
ARLINGTON,  VIRGINIA  22202 

OR.  E.  BAUER 


UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA.  15213 

OR.  N.  2A0USKY 
DR.  M.  BIOnOI 
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